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Abstract 



The quantum lightcone fluctuations in flat spacetimes with compactified spa- 
tial dimensions or with boundaries are examined. The discussion is based 
upon a model in which the source of the underlying metric fluctuations is 
taken to be quantized linear perturbations of the gravitational field. General 
expressions are derived, in the transverse trace-free gauge, for the summation 
of graviton polarization tensors, and for vacuum graviton two-point functions. 
Because of the fluctuating light cone, the flight time of photons between a 
source and a detector may be either longer or shorter than the light prop- 
agation time in the background classical spacetime. We calculate the mean 
deviations from the classical propagation time of photons due to the changes 
in the topology of the flat spacetime. These deviations are in general larger in 
the directions in which topology changes occur and are typically of the order 
of the Planck time, but they can get larger as the travel distance increases. 
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I. INTRODUCTION 



The existence of fixed lightcone structures is one of the characteristics of classical gravi- 
tational theory. Lightcones are basically hypersurfaces which distinguish timelike separation 
from spacelike separation and divide spacetime into causally distinct regions. However, if 
gravity is to be quantized, it is natural to expect that the quantum metric fluctuations would 
smear out the lightcone, and the concept of a fixed lightcone structure has to be abandoned. 
Based upon the observation that the ultraviolet divergences of quantum field theory arise 
from the light cone singularities of two-point functions, and that quantum fluctuations of 
the spacetime metric ought to smear out the light cone, thus possibly removing these sin- 
gularities, Pauli im conjectured many years ago that the ultraviolet divergences of quantum 
field theory might be removed if gravity is quantized. This idea was further explored by 
several other authors P-^. At present time, this conjecture remains unproven. If light- 
cones fluctuate, so do horizons, which are, of course, lightcones. The horizon fluctuations 
could then presumably lead to information leakage across the black hole in a way that is 
not allowed by classical physics. Bekenstein and Mukhanov have suggested that horizon 
fluctuations could result in discreteness of the spectrum of black holes. Since the existence 
of black hole horizons is the origin of the so-called black hole information paradox, which 
has been widely discussed in the literature but still remains to be resolved, the study of 
lightcone fluctuations might help us better understand the problem. 

Recently the problem of light cone fluctuations has been investigated in a model 
of quantum linearized theory of gravity, where the fluctuations are produced by gravitons 
propagating on a background spacetime. The lightcone is smeared out if the linearized 
gravitational perturbations are quantized. It has been demonstrated that gravitons in a 
quantum state, such as a squeezed vacuum state, or a thermal state, can produce light 
cone fluctuations, thus smearing out the light cone. Because of the fluctuating light cone, 
the propagation time of a classical light pulse over distance r is no longer precisely r, but 
undergoes fluctuations around a mean value of r. The fluctuations in the photon arrival 
time can also be understood as fluctuations in the velocity of light. This model has been 
applied to study the quantum cosmological and black hole horizon fluctuations [||. It is 
interesting to note that recently, the quantum gravitational metric fluctuations have also 
been discussed within a different context, i.e., a Liouville string formulation of quantum 
gravity [^|10| . In this paper we shall examine light cone fluctuations in flat spacetime with 
nontrivial topology based upon the model proposed in Ref . In Sec. II, we review the 
basic formalism and examine its gauge invariance, then derive general expressions for the 
vacuum graviton two-point functions in the transverse trace-free gauge. In Sec. Ill we study 
the light cone fluctuations in flat spacetimes with a compactified spatial dimension, and with 
a single plane boundary. Our results are summarized and discussed in Sec. VI. 

II. BASIC FORMALISM AND GRAVITON TWO-POINT FUNCTION IN 
TRANSVERSE TRACE-FREE GAUGE 

Let us consider a flat background spacetime with a linearized perturbation hfj,u propa- 
gating upon it , so the spacetime metric may be written as 
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ds^ = g^udx^dx" = (r/^j, + h^i,)dx^dx^ = dt^ — dy^ + h^^dx'^dx" . (1) 

Let a{x, x') be one half of the squared geodesic separation for any pair of spacetime points x 
and x', and <Jo{x, x') be the corresponding quantity in the flat background . We can expand, 
in the presence of the perturbation, a{x, x') in powers of h^^, as 

a = (To + o"! + o"2 H , (2) 

where cxi is first order in h^^, etc. We now suppose that the linearized perturbation h^j.^ is 
quantized, and that the quantum state is a "vacuum" state in the sense that we can 
decompose h^j^y into positive and negative frequency parts h^^ and h~y, respectively, such 
that 



/i+j^) = o, (^|/^;. = o. (3) 

It follows immediately that 

(V) = (4) 

in state \ip). In general, however, ((/i^zy)^)/?, 7^ 0, where the expectation value is understood 
to be suitably renormalized. This reflects the quantum metric fluctuations. 



A. Basic formalism and gauge invariance 



If we average the retarded Green's function, Gret(x, x'), for a massless scalar field, over 
quantized metric fluctuations, we get 

(^-^^'^')) = ^^/^ '"P(-2R)) • 

This form is valid for the case in which (cr^) > 0. It reveals that the delta-function behavior of 
the classical Green's function, Grei, has been smeared out into a Gaussian function peaked 
around the classical lightcone. This smearing can be understood as due to the fact that 
photons may be either slowed down or speeded up by the light cone fluctuations. Photon 
propagation now becomes a statistical phenomenon, with some photons traveling slower 
than the light on the classical spacetime, and others traveling faster. Note that the Gaussian 
function in Eq. (|^) is symmetrical about the classical light cone, do = 0, and so the quantum 
fluctuations are equally hkely to produce a time advance as a time delay. 

Light cone fluctuations are in principle observable. It has been shown, by considering 
light pulses between a source and a detector separated by a distance r, that the mean 
deviation from the classical propagation time is related to (a^) by P| 




Note, however, that At is the ensemble averaged deviation, not necessarily the expected 
variation in flight time of two photons emitted close together in time. The latter can be much 
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smaller than At due to the fact that the gravitational field may not fluctuate significantly in 
the interval between the two photons. This point is discussed in detail in Ref. 0. In order 
to find At in a particular situation, we need to calculate the quantum expectation value 
(af) in a chosen quantum state. For this purpose, we first have to compute cxi for a given 
classical perturbation along a certain geodesic, then average over the quantized metric 
perturbation. If we consider a null geodesic specified by 

df = dx^ - h^^dz'^dx", (7) 

then by following the same steps as those of Ref. 0, we can show that in a general gauge 

ai = lAr r h^^n^'n'^dr, (8) 

Z Jro 

and 

{crD = hArf r dr T dr' n^^n'^n^n'^ {KAx)h,„{x')) r . (9) 

Here dr = \dx\, Ar = ri — Tq and n'^ = dx^ /dr. The graviton two-point function, 
{hf_iu{x)hp^{x'))R, is understood to be renormalized, so that it is finite when x = x' and 
vanishes when the quantum state of the gravitons is the Minkowski vacuum state. 

A few comments on the derivation of Eq. (|^) are in order here. It is obtained by averaging 
the Fourier representation of a 5-function. It may come as a surprize that although we started 
with an analytic expansion of a in powers of h^^,, the result is not analytic as (a^) 0. 
This arises because we use the first order expansion of a in the argument of an exponential 
function, but afterwards retain all powers of h^^. One can reasonably ask whether this is 
a valid procedure. A test of the self-consistency is to retain the (T2 term and then follow 
the same procedure. The result is Eq. @ with ctq replaced by cJq + o^- This has the same 
physical interpretation as before; the only effect of the 02 part is to shift the location of 
the mean lightcone. Thus in this order we encounter the backreaction of the gravitons in 
perturbing the original classical geometry to a new classical geometry. Although this is less 
than a complete demonstration of the validity of Eq. (^, it does indicate that it arises from a 
self-consistent calculation. In any case, the only result that we really need in the remainder 
of this paper is Eq. (^, which may be derived either from Eq. (^), or else more directly by 
averaging the square of Eq. (|) . 

Let us now turn to the question of the gauge invariance of the formalism. Under a gauge 
transformation specified by 

x'^ = x^ ^i^[x), (10) 

^i,v{A = h^u{x) - ^m{x) , (11) 

where C,^{x) is of order hfj_^, the quantities ui and (af) are not in general invariant. However, 
we can show that this is due to the fact that At is a coordinate time interval rather than a 
proper time interval. To better understand the gauge invariance, let us examine a situation 
in which a light signal travels between two points in space labeled by P and Q, with a 
classical metric perturbation h^i, in the intervening region, as illustrated in Fig. 1. 
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FIG. 1. A light ray ( dashed hne ) makes a round trip travel between two points, P and Q, in 
space. 

For simplicity, let us assume that the propagation is in the x-direction. We shall look at 
the travel time in two different gauges, or coordinate systems, primed and unprimed. For 
light rays traveling in x direction, we have 



dx V dx dx 



(12) 



Here the upper sign is used for outgoing light rays and the lower sign for incoming rays. So, 
one way travel time St in the unprimed gauge is 



St 



P-*Q 



dx^ dx" , 
X ] — — - — dx 



j dx - \ j V( 

Jxp Z Jxp 

fXQ 1 r« , ' / ,x dx^' dx" , 1 r« ^ /X dxf^ dx" , 



(13) 



which, within the linearized theory, can approximated as 



/■^o , 1 /-^Q , ..dx'^dx" , , I ^ , .dx^dx" , 
St,^^ = dx - - J dx - - dx, 



rXQ 1 ri 

Jxp Jx, 



. 1 /"^Q , / , ,, dx^" dx" 
A J X p dx dx 



dx' — / dx — I dx 



^« d^ 

xp dx 



xp 



dx 



XQ{t') - UQ, t') - ( ^p{to) - UP, to) ) - CtiQ, t') + it{P, to) 
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-- r h' (x')^^dx' 

2 7.; ''''^^^^dx' dx' "^"^ ' 

(14) 

where we have used the fact dt/dx — 1 for outgoing Hght rays within our approximation. 
Similarly, we have 

,xp I , , . dx'' dx'' , 

dtg^P = - / dx +- h^„{x) — — dx 
Jxq ^ '^^Q ^x dx 

1 f^Q , , dx'' dx'' 1 f^p dx'' dx" 



r^Q 1 f-Q , dx''dx'' , f-p di^ f-p dit , 
/ dx — - I n [x j—TY—rj- dx + -r— dx — — dx 
Jxp ^ 'xp dx cLx dx ^ dx 

XQ{t') - UQ, t') - ( xp{Q - UP, Q ) + t') - Ct{P, Q 

-- r'^h' (x')^^dx' 

2 7.' "^"^""Ux' dx'"^""' 



(15) 



using the fact that for incoming light rays, dt/dx — —1. Note that 

x'Q{t)^XQ{t)-UQ,t), (16) 



x'p{t)^xp{t)-UP,t), (17) 

so 

5tp^Q = St'p_^Q - CtiQ, t') + 6(P, to) , (18) 

and 

5tQ^P = St'Q^p + it{Q. t') - it{P. t'o) . (19) 
It follows that the round trip travel time is 

At - 5tp^Q + 5tQ^p = At' + et(P, to) - UP, t'o) ■ (20) 

Therefore, the one way travel times, Stp^g and Stg^p are, in general, not invariant unless 
both the source and the detector are outside the regions where gravitational perturbations 
hij,,y are non-zero. In that case, it is physically reasonable to set ^{P,t) and ^{Q,t) to zero. 
Similarly, the round trip time At is invariant only if the source ( it also acts as a detector 
in this case ) is outside of the gravitational perturbations. 

However, it is interesting to note that the round trip proper time interval for the source. 
At, is gauge invariant. Denote the proper time intervals in two different gauges by At and 
At', and keep in mind the fact that on the world line of the source, generally, ^ << 1. We 
then have 
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At' = J ^/l + h^odt' = Jdt' + ^J h'oodt' 

= At' + to) - UP, Q + f^oodt 
= At + ^J hoodt = AT, 

(21) 

where we have used Eq. (pOl). This shows that we should really consider how proper time 
rather than the coordinate time is affected by light cone ffuctuations. However, the calcu- 
lation of the proper time in a general gauge is a rather difficult task, because the source 
(and detector) may not be at rest with respect to the chosen coordinate system, and thus 
in general the emission and the subsequent reception may not happen at the same point 
in space. To find the proper time, we have to integrate along the geodesic between two 
events, the emission and the subsequent reception. In general, there is a Doppler shift due 
to fluctuations in the positions of the source and the mirror. However, the analysis can be 
greatly simplified if we adopt the transverse-tracefree ( TT ) gauge, which is specified by 
the conditions 

= djh'^ = h^'' = 0. (22) 
To see this, let us examine the geodesic equations for a test particle 

d'^x'^ dxP dxP . 

which, when written in term of derivatives with respect to coordinate time t, becomes 

d'^x^^ ^„ dxP dx'^ ^. dx^ dxP dx" „ 

dH dt dt dt dt dt ' ^ ^ 

For a non-relativistic test particle, ^ << 1, so, to the leading order. 



dt 



dH 



■ (25) 



But, in the TT gauge, FJ^ = 0. Therefore, from the above equation, we can see that 
if the test particle is at rest at t = 0, then it will subsequently always remain at rest |]13 



So, if we are considering the emission and reflection of a light signal between two points 
(particles) in the TT gauge, then the proper time 6t between emission and reception (after 
reflection) of the signal is related with the coordinate time by 

= J VSttdt = J ^(1 + hoo)dt = J dt = 6t. (26) 



Here we have appealed to the fact that /loo = in the TT gauge. Therefore, the coordinate 
time for the round trip in the TT gauge is the proper time, and At calculated from Eq. (^ 



7 



in the TT guage is actually a gauge invariant quantity. In this gauge, the mean squared 
fluctuation in the geodesic interval function reduces to 

{al)= ]{ArY Hdr H dr' n'n^n^rT {K,{x)hkM))R 

= \{Arf r dr r dr' n'n^n^rT {h,,{x)hm{x') + h,,{x')hm{x))R • (27) 

Here = dx^ /dr is the unit three- vector defining the spatial direction of the geodesic. 



B. Graviton two-point function in transverse trace-free gauge 

If we work in the TT gauge, the gravitational perturbations have only spatial components 
hij and they may be quantized using a plane wave expansion as 

= E KAe.,(k, A)/k + H.c]. (28) 

k,A 

Here H.c. denotes the Hermitian conjugate, A labels the polarization states, and 

/k = (2cu(27r)3)-5e^(k— -*) (29) 

is the mode function, where 

w=|k|, \k\ = {kl + kl + klY^, (30) 

and the e^,y(k. A) are polarization tensors. ( Units in which 327rG' = 1, where G is Newton's 
constant and in which h = c = 1 will be used in this paper.) 

Now we shall first calculate the Minkowski spacetime Hadamard function for gravitons 
in the transverse tracefree gauge. Let us define 

G%i^,x') = {0\h,,{x)hkiix') + h,^{x')hMixm . (31) 

Then we have 

G^ii^, = 1^ / rf'kX: e.,(k, A)e,Kk, A)i-e^'^ (— ')e-(*-*') . (32) 

Equation ([A8|) in the Appendix for the summation of polarization tensors in the transverse 
tracefree gauge gives 

Cijik, \)eki(k, A) = SikSji + SuSjk - SijSki + kikjkkh + kikjdki 

X 

-\-kkki6ij kikiSjk kik^Sji kjki6ik kjk^Su , (^3) 

where 

k = ^. (34) 
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We find that G% {x, x') can be expressed as |[Tl| 



G%ix,x')= 2Re i5,kS,i + 5u5,k - S,,5ki + A,) x -^^ j rf3i^i_g^k.(x-x'),-.c.(t-t') 



e 



(2vr) 

2Re {5ikSji + - 6ij6ki + Ai) x (O|0(x)0(x')|O) , (35) 



where we have defined a formal operator 



'd^d' d,di a.a' 9.9' 9,9' 9,9' a9^^ 



and (O|0(x)0(x')|O) is the usual scalar field two-point function. Here the formal operator V^^ 
should be understood in the sense of a Green's function, but when we do our calculations 
in momentum space its effect is to bring in a factor of fc"^. 
The combination of these results with Eq. (|2^ ) gives 



{al) = -{ArfJ^ drj^ dr' 1-^ >- + ^ ^ (0(a;)0(x'))R. 



Introduce two functions Fij{x,x') and Hijki{x,x') by 

V2 



Fij{x,x')= i?e^(O|0(a;)0(a;')|O) 



and 



9- 9' 9 9' 

A,m(x,x')= i?e ' {0\(Pix)<P{x')\0) 

mOkd' 1 



1 . 3k^e^k.(x-x')^-..(.-0 



(27r)3 J 2uj^ 



gI^-Ii can be expressed as 



^ijki — '^Fij^ki + 2FkiSij — 2FikSji — 2Fii6jk — 2Fji6ik — 2Fjk6ii + 2Hijki 



(37) 



V2 (27r)3i 2CU 
(27r)3y ''^2a;3^ ' ' ^^^^ 



V4 (27r 

Re f ^3^^ fcjfc,fcfcfc; ^ik-(x-x')^-»a;(t-t') _ ^3g^ 



+2^(1) (x, x')(5ifc(5,7 + - (40) 

where 



is the usual Hadamard function for massless scalar fields with 2(7 q = (t — t')' — (x — x' 



t'f _ (x - ^'^2 

and Fij{x,x') and Hijki{x,x'), which will be calculated in the Appendix, are given by 

2" 



(2vr) 



li ...2 .2, R + At 

- \n(R^ - At^) + — In — 

2 ^ ^ 4R \R-At 



(42) 



and 



Hijki{x,x') = T^^did'jdkd'i 



{R^ + SAt^) ln(i?=^ - At 



2\2 



+ (sRAt 



At^' 
R 



In 



R + At 
R- At 



(43) 



Here R= |x — x'l. 



III. LIGHTCONE FLUCTUATIONS IN FLAT SPACETIME WITH NONTRIVIAL 

TOPOLOGIES OR BOUNDARIES 

In this section, we study lightcone fluctuations in two cases: flat spacetime with a com- 
pactified spatial section, and with a single plane boundary. 



A. Flat spacetime with a compactified spatial section 

Let us now assume that the spacetime is flat but compactified in the z direction with 
a periodicity length L ( "circumference of the universe" ). This means the spatial points z 
and z + L are identified. The effect of the space closure is to restrict the field modes to a 
discrete set 

/k = (2Lu(27r)2L)-5e^(''---*) (44) 

with 

kz = ^, n = 0,±l,±2,±3,... (45) 

We now analyze the lightcone fluctuations, assuming that the gravitons are in the new 
vacuum state \0l) associated the discrete modes of Eq. ). 

First consider a light ray along z direction, i.e. along the direction of compactification 
(Fig. 2 ), propagating from point (0,0, a) to point (0,0,6) in space, then, we have from 
Eq. (P 

{(^!)=lib-af f dz f dz' {QL\Kz{x)hzz{x') + hzz{x')hzz{x)\QL)R 
o J a J a 

= \(h-af fdz f'dz'G^^l1it,0,0,z,t',0,0,z'). (46) 

O J a J a 
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Here we have defined 

= {QL\hUx)h,,{x') + K,{x')K,{x)\Ql) - {Q\Kz{x)K,{x') + h,,{x')h,,{x)\Q) , (47) 
and the integral is to be carried out along the geodesic. 



(t,0,0,a) 



(t ' ,0,0,b) 



FIG. 2. A light ray ( dashed Une ) propagates in the direction of compactification in a 
cylindrical "universe" from point (t,0, 0,a) to point (t', 0,0,6). Here only two spatial dimensions 
are plotted. 

If we adopt the notation 

(t,0,0,z, t',0,0,/) = (i,^, i',^'), (48) 
the renormalized two-point function can be found by using the method of images to be 

+ 00 ' 

G(i^(t,z,t',/)= E GiL(t,^,t',^' + nL) 

n=—oo 

+ 00 ' / 

= 2 ^ {D^^\t,z,t\z' + nL)-2F,,{t,z,t\z' + nL) 

n=—oo \ 

(49) 

where the prime on the summation indicates that the n = term is omitted. Substituting 
Rt = into Eq. ( |A25| ) in the Appendix and replacing Ax by Az, we have 



n +00 



n=— oo 

2 



+ 



At' 



In 



{Az-nLY 4{Az-nL)^ \Az -nL + At 



Az-nL- At 



At 



In 



3{Az-nLy A{Az-nL)^ \Az-nL + At 



Az-nL- At 



(50) 
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For the null geodesic 



At = Az, 



(51) 



we get, after an evaluation of the integral, 



dz / dz'Gg\1{t,Z, t',z')\^t=Az 
J a 

k^{n^-2e^) (n + 2e)^ 
^2)^^ 2(n + e)3 

{n + 2ef 



E 

n=—OD 
+00 



16e' - 2e' 



127r2 
1 

12^ 
1 



n 



2 n 



n 



E 

n=l 

+00 

E/(^'e) 

n=l 



(n^ — e2)2 



+ 



in + e) 



n 



where we have defined 



(b-a) 



r 



and 



fn^ — e2)2 



+ 



n 



n 



(52) 



(53) 



(54) 



It appears that there is a singularity in the summand f{n,e) whenever n = e, i.e., 
whenever the distance r is an integer multiple of L. However this singularity is illusionary, 
as it should be from a physical point of view since there is nothing special when n = e. This 
can be seen if we expand the summand at the point e = n to get 



X/ N 19 27, 271n(3) + 68, , 

f{n, e)^- + - In 3 + Y ^ " ^ + ^ ^ " ^ • 

3 4 8n 



(55) 



So, f{n,e) is finite as e approaches n. Note also that 2e = n is also not a singularity. The 
summation converges, as the asymptotic form of f{n.t) as ri ^ 00 is 



32e2 



+ 0{n 



(56) 



However, a generic closed form result for the summation is hard to find. So we now discuss 
two special cases. The first is the one in which the distance traversed by the light ray is 
much less than the periodicity length, h — a <^ L. Then we get 



j-b i-b 

/ dz / dz'GTAx,x')^Y: 
Substitution of this result into Eq. (H3) yields 



^2 1 



4e^ 

37r2 n2 9 



(57) 
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Therefore the mean deviation from the classical propagation time is 



Since we are working in Natural Units, this result reveals that the mean deviation in 
travel time is less than the Planck time and grows linearly with increasing r when r is small 
compared to the periodicity length L of the universe. 

If e 1, i.e., r ^ L, the light loops around the " universe ", and summation Eq (|52|) 
can be approximated by the following integral 

16{x^ - 2 " 



+ 



(x^ - 11 



(60) 



Evaluating the integral with the aid of the computer algebra package Maple, series expanding 
the result and keeping the leading terms only, we arrive at 

81n(2e 

Therefore the mean deviation from the classical propagation time is 



dz dz'GTAt,z,t',z')^e-"—^. (61) 



w (c^i ) 1 IT 

where r is assumed to be much greater than L. So the lightcone fluctuations can, in principle, 
get as large as one would like if the light ray travels around and around. This is interesting 
in the sense that it suggests that a fluctuation which is much greater than the Planck scale 
could be achieved. 

Now we turn to the case where the light ray moves along the direction perpendicular to 
that of compactification, for instance, along x direction. If the light ray travels from point 
(a, 0,0) to point (6,0,0), as illustrated in Fig. 3, then 

1 fb fb 

{(^i)=zip-af dx dx' {OL\h^^{x)h^^{x') + h^^{x')h.,^{x)\OL)R 



l{b-af [''dx dx'G^^^g{t,x,0,0, t',x',0,0), 

J a J a 

1 pb nb ' 

^{h-af dx dx' G^^L{t,x,0,0,t',x',0,nL). (63) 

o J a J a 
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i (t,a,0,0) 



FIG. 3. A light ray ( dashed hne ) propagates perpendicular to the direction of compactifi- 
cation in a cylindrical "universe" from point (t,a, 0,0) to point (t',b, 0,0). Here only two spatial 
dimensions are plotted. 



Let us now define 



p = X — x', b — a = r, 



(64) 



then if we use Eq. ( |A25D in the Appendix and bear in mind the fact that for the hght ray 
At = Ax, we have 



Gii'Kt, X, 0, 0, t', x', 0, nL) = g,{p) + g^ip) 



(65) 



where 



91=2 j: 



1 p2 {riLY 



+ 



47 p^ {nLf 



87r2 (p2 + (nL)2)4 37r2 (p2 + (nL)2)4 127r2 (p2 + (nL)2)4 



(66) 



and 



+ 00 ' Y 

92 = -2 y — ^ In 



p2 + (nL)2 + p 



^^p^ + {nLf-p^ 



p{nLf 



p3 [riLY 



16 (p2 + (nL)2)(9/2) 4 (p2 + (nL)2)(9/2) 



p5(nL)2 



2 (p2 + (nL)2)(9/2) 



We can clearly see is an even function of p, so, 



fdx rrfx'GWi?(t,a;,0,0, t',a:',0,0) = 2 f dp{r - p){g, + g^) . 

J a J a JO 



(67) 



(68) 



Performing the integration (integrate by parts for those terms involving logarithmic func- 
tion), we arrive at 
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2 / dp{r-p){gi + g2) 
Jo 

^ ' r 



9 +00 



2(e2 + n2)2 4(e2 + n2)2 



24(n2 + e2)5/2 



In 



\/n? + e'^ — e 



(69) 



where e = r/L as before. The above series can be shown to be convergent. Yet a result in 
closed form is not easy to find. Let us first examine the case in which r <^ L, where 



pb rb 

/ dx / cia;'G«f,(x,x')^E 



457r2 452 x 21 



Here we have used 



+00 1 



TT 



n=l 



45 X 21 

Thus the mean deviation from the classical propagation time is 

2 27r2 



(70) 



(71) 



At 



21 45 VL 



(72) 



This result holds in the small e regime. The time deviation is much smaller than that 
for light rays propagating along the compactification direction ( compare with Eq. (^) ). 
This reveals that light cone fluctuations due to topology change are more likely to be felt 
in the direction of compactification than in the transverse direction, if we perform local 
experiments in which r, the distance between the source and the detector, is very small as 
compared to L, the periodicity length. 

We now turn our attention to the case in which r ^ L, i.e., e ^ 1. Here it is easy to see 
that the summation in Eq. (B^) can be approximated by the following integral 



f' dx f dx' G'-^l^.it, X, 0, 0, t\ x', 0, 0) 

a J a 



4e /■°° , 
/ dx 

Jl/e 



X 



+ 



2(1 + X2)2 4(1 + X2)2 



In 



24(x2 + 1)5/2 \^^/^2^_l^ 



(73) 



If we perform the integral and series expand the result, we have, to the order of 0(e) 

r dx f dx' G^^l^^it, X, 0, 0, t', x', 0, 0) 

J a J a 

^ cje - 4 ln(e) , 



(74) 
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where ci and C2 are constants given, respectively, by 



:? = / dx 

Jo 



ln(x + Vx^^TT) 



2.468 . 



(75) 



X\^X^ + 1 



and 




(76) 



and the mean deviation from the classical propagation time i 




(77) 



_ ci 1 rr 



(78) 



This result applies in the regime where r ^ L. Here we have the same functional dependence 
on r as in the case where light rays loop around the compactified dimension many times. The 
only difference lies in the proportionality constants. In fact, here the mean time deviation 
is also smaller than that for light rays traveling in the direction of compactification, since 
the numerical constant Ci/vr ^ 0.5. 



Let us assume that there is a single plane boundary located at z = in space such that 
metric perturbations satisfy the following Neumann boundary condition (The reason that 
we use the Neumann boundary condition instead of the Dirichlet boundary condition here 
is to get a positive (erf). ) 



where kt and xt denote the components of k and x, respectively, in directions parallel to the 
boundary. Now if we assume that the gravitons are in the vacuum state |0') associated with 
the modes of Eq. (pO|), we have, for a light ray propagating perpendicular to the boundary 
from point (0, 0, a) to (0, 0, b), 



B. Single plane boundary 





(80) 



{al)=-{h-af dz dz' {Q'\h,,{x)K,{x') + K,{x')K,{x)\Q') 

O J a J a 



1 rb rb 



R 




(81) 
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(t,0,0,a) 



(t ,0,0,b) 



FIG. 4. A light ray ( dashed hne )propagates in the direction perpendicular to the plane 
boundary, starting a distance away from the boundary 

Here the renormalized graviton two point function G^J-]^{x,x') can be found by the 
method of images as usual and the only difference is an overall sign change as we go from 
the Dirichlet boundary condition to the Neumann boundary condition. The reason for this 
is that, to satisfy the Neumann boundary condition, we need to add the image term instead 
of subtracting it as in the case of the Dirichlet boundary condition. So, G^^J^{x,x') may be 
obtained by picking out the n = term in Eq. (|^) and setting Az = z + z' to get 



2(t - 1 



!\2 



it - t') 



n'^iz + z'Y 
4 



+ 



2Tx'^{z + z'f 

it - t') 



In 



In 



' z + z' -{t-t') \ 
^z + z+{t + t')) 

' z + z' - [t-t'Y 



37r2(z + ^')2 27r2(;z + z')3 + z' + {t + t') ) 

Substituting this result into Eq. ( |STD and performing the integration, we finally get 

{h - af [62 _ a2 + (a2 + Aah + fe^) \^^k] 



(82) 



247r2(6 + a) 



(83) 



Note that this result is always greater than zero. However, had we chosen the Dirichlet 
boundary condition, we would have that (cr^) < 0. Recall that the formalism which we are 
using applies only if {af) > 0. 

When the hght ray starts very close to the boundary such that a <^ r, we have 



The mean deviation in travel time is 

At = 



247r2 



;i + ln(r/a)) . 



1 + ln(r/a) 



(84) 



(85) 



r V 247r2 

which diverges as a approaches 0. This is not surprising since the energy density of a 



quantized field blows up on the boundary. However, it has been shown recently Il2| that, 
if one treats the boundaries as quantum objects with a nonzero position uncertainty, the 
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singularity in energy density is removed. The result, Eq. (|85D, applies whenever r 3> a. 
The other limit is when r <^ a, where the mean squared fluctuation in the geodesic interval 
function is approximated as 



^4 



consequently, the mean deviation in time is given by 



(^i) ^ TTTT^ , (86) 



At = ^^ = ^^-. (87) 

We now consider a null geodesic which is z distance away from and parallel to the plane 
boundary. The relevant renormalized Hadamard function is given by Eq. ( |A25D with !S.z 
being replaced by 2; + 2;'. 



A (t ,b,0,z) 



1 (t,a,0,z) 



FIG. 5. A light ray ( dashed line )propagates in the direction parallel to the plane boundary, 
starting z distance away from it 



Now suppose the geodesic starts at point (t,a,0,z) and ends at point 0,2;), then 
the mean squared fluctuation in the geodesic interval function is 



(^?) = \{^-(^f tdz fdz' Gf,,,(t,a;,0,z, t',a;',0,z). 



a J a 



(88) 



Here x, 0, 2;, t',x',0,2;) is also given by Eqs. ( |65| - |67|) but with a replacement of riL 

hj 2z. Therefore 



b j-b 

dx / dx' G)^l^^{t, X, 0, z, t\ x\ 0, z) 

a J a 

2 



xxxx \ 
,4 



2(e2 + 4)2 (e2 + 4)2 



8e5 + 32e3 + 48e , f + ^ + e 
H :rTT-. . In ' 



24(4 + e2)5/2 \^74^_g^ 



(89) 
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where e = rjz. Since the above expression is very comphcated, we shall discuss two inter- 
esting special cases. One is when r ^ then we have 



{a')^—ln{r/z) (90) 



and 



V^_/M^ (91) 

This also blows up as z approaches 0, however the functional dependence upon z is different 
from that of Eq. (B5|). The other limit is when r -C 2;. For this case, we find 



(a^) ^ (92) 



and 



At = xi^ = 1 . (93) 



12v^7r \z 



IV. SUMMARY AND DISCUSSION 

In this paper, we have obtained general expressions, in the transverse tracefree gauge, for 
the vacuum graviton two-point function for various boundary conditions. These were used 
to study the lightcone fluctuations in fiat spacetimes with a compactified spatial section and 
with a plane boundary. The mean squared fluctuations of the geodesic interval function and 
therefore the mean deviations from the classical propagation time have been obtained. 

In the case of a compactified spatial section, when the travel distance is less than the 
periodicity length, the fluctuation in the propagation time is less than the Planck time. 
In this limit, the effect is much larger for propagation in the periodicity direction than 
for propagation in the transverse direction. Thus the local lightcone fluctuations become 
anisotropic, reflecting the global structure of the spacetime. When the travel distance is large 
compared to the periodicity length, the fluctuation in travel time increases with the square 
root of the distance traveled for propagation in either direction, and the only difference lies 
in the proportionality constants. Here we have a possibility of having fluctuations larger 
than Planck scale by several orders of magnitude. 

In the case of a plane boundary, as light rays start closer and closer to the boundary, the 
lightcone fluctuations blow up as the square root of the logarithm of the starting distance 
both when light rays propagate perpendicularly and parallelly to the boundary. This is 
not as surprising as it might seem because the imposition of a fixed boundary can lead 
to singular expectation values of local observables, such as energy densities. However we 
expect this singularity to disappear if one treats the boundary as a quantum mechanical 



object with a nonzero position uncertainty |T^. It is also found that if the starting distance 



from the boundary is fixed, then the fluctuation in travel time grows as the square root of 
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the logarithm of the distance traversed when this distance is large compared to the starting 
distance. 

In summary, we have demonstrated that in the linearized theory of quantum gravity, 
changes in the topology of flat spacetime produce lightcone fluctuations. These fluctuations 
are in general larger in the directions in which topology changes occur and are typically of 
the order of Planck scale, but they can get larger for path lengths large compared to the 
compactification scale. It is interesting to note that this effect could become significant in 
theories which postulate extra dimensions compactified on a very small scale. 
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APPENDIX 



A. Summation of graviton polarization tensors in the TT gauge 

Let us introduce a triad of orthonormal vectors (ei(k), e2(k), e3(k)) with 

e3(k) = ^ = k, (Al) 

the unit vector in the direction of propagation. The triad satisfies the orthonormahty relation 

ea{k) ■ eb{k) = 6ab, a, 6 = 1,2, 3. (A2) 

This relation can be written, in terms of the components in the coordinate system charac- 
terizing the metric, as 

el{k)el{k) = 6ab, a, 6 =1,2,3. (A3) 

Here the Einstein summation convention is employed. We also have 

ej.(k)e^,(k) = e\e{ + e^ei + k'^ = i, j = x, y.z. (A4) 

Therefore, the two independent graviton polarization tensors in the TT gauge are given, in 
terms of the triad, by 

e*^(k, +) = ei(k) ® ei(k) - eUk) ® ei{k) , (A5) 
e'' (k, X ) = e\ (k) ® ei (k) + (k) ® e{{k) , (A6) 

(A7) 

where we have adopted the notation of Ref [|13|. Hence, 

J2 eij{k, A)efc/(k, A)= Cijik, +)eki{k, +) + eij{k, x)efc/(k, x) 

A 

id''' ID^ tD^ tJ' /D^ /D^ /D^ t:>^ _L /^-^ ^ /^^ 

— C^CiCidi ^'1^2^\^\ ^ ^2^2^2^2 

p}^ J- p^ fD^ P^ P^ -1- P^ P^ fD^ I .pfe 

= (eiet + 6*262 )(ei el + 6^362) + (eie^ + e^e'2)(eie^ + 6^262) 

-(elei + e^e^2)(eM + 44) 
= {5'^ - k'k^){5^^ - + (5^' - k'y){6^'' - Vk^) 

= ^ik^jl + ^il^jk — ^ij^kl + kikjkkki + kikjSkl + 

kkhSij - kikSjk - kihSji - kjki5ik - kjkk^u ■ (A8) 
This result can also be obtained as follows. Let us introduce a 4th-rank tensor 

T'^^^{k) = e''{k, A)e'='(k, A) , (A9) 

A 

which has the following symmetry properties 



21 



rpijkl rpjikl rpijlk rpklij ^AlO) 

However, the objects, which are at our disposal to construct T^i^^^ are only /c* and 5^^ , thus 
in general, we have 

(All) 

where A, B, C, D, E are constants to be determined. This tensor is subject to the transver- 
sality condition 

and the trace-free condition 

rpiikl _ rpijkk _ g (A13) 

Applying these constraint conditions to T'^^^'- and solving the resulting equations leads to 

a = d = -e = -c = -b . (A14) 

Therefore T*-^'^' is the same as the right-hand side of Eq. (XH), apart for a multiplicative 
normalization constant which can be chosen to be unity. 



B. Vacuum graviton Hadamard function in the TT gauge 

Here we evaluate the function Fij{x,x') and Hijki{x,x') defined in Eqs. (^) and (^3D, 
respectively. Once these functions are given, the graviton two point functions are easy to 
obtain. Define 



R= ^[x-x'y + {y-y'y + {z-z'y, At = t-t', k = \k\ = u . (A15) 

Then, 



J?..(r r'l- [ ^3i^^ ik.(x-x')„-ia;(t-t') 

t,,[x,x)- J (^^^^^e e 



(27r 

T}p foo p—ikAt ^71- „27r 



(27r)3 ' ^ Jo 2k Jo Jo 

I ^ ^. I dk 

:OiO, — I TT7 sin kR cos A; At . 



(A16) 



Because there is an infrared divergence in the above integral, we will introduce a regulator 
(5 in the denominator of the integrand and then let (3 approach after the integration is 
performed. 
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2 2 r+oo f^j^ 

Fij{x,x')= . .^ did', lim — / sin kR coskAt 

(27r)2 ^ f3~^o RJo k^ + 



p^aS^lim/(/3.fl.A<). 



(A17) 



Here we have used a integral in Ref. |n| and defined 

R, At)= -i-|e^(^*-^)Ei[/5(i? - At)] + e-^(^*+^)Ei[/?(i? + At)] 
AfjR [ 

-e^(^*+^)Ei[-/3(/? + At)] - e^(^-^*)Ei[/5(At - i?)]| . 



(A18) 



Here Ei(x) is the exponential- integral function. Making use of the fact that, when x is small, 

(A19) 

where 7 is the Euler constant, and expanding / around /3 = to the order of we get 



1 1 

Ei(x) ^ 7 + In Ixl + X + -x^ H x^ + Oix'^) 

4 18 



/(/3,i?,At)^l-7-ln/3--ln(i?2-At2)-^ln''^ + ^^» ' ^"^-"^ 



m \R- At 



(A20) 



Taking the limit and keeping in mind that the constant terms ( with respect to x and x' ) 
vanish under differentiation, we finally obtain 



Fij (x, X ) 



1 



(2vr; 



1 , / . 2^ At fR+ At 

- In - At^ + — In — 

2 ^ ' AR \R-At 



(A21) 



Now let us turn our attention to Hij^i^x^x'). We have, proceeding with similar steps as 
we did for Fij{x, x'), 



Hijkiix, x') 



Re 
(2^ 
1 



r ,3 fci fcj fcfc fc; ik- (x-x' ) - Af 
J 2^5 



:did'.dkd'i lim 



1 d r+ 



dk 



2I3R d(3 Jo fc2 + /32 



(27r)2 ' ^ 

(^Wlimi^|miS,A*). 



■ sin kR cos A; At 



(A22) 



Now expand jg^fiP, R, At) to order to find 



2/3 a/3" 



(ln/3 + 7 - l)i?^ + 3(ln /3 + 7 - 1) At^ 



12R 



{R + At)^ In \R + At| + (i? - At)^ In \R - At\ 



(A23) 
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Plugging this result into Eq. (A22) and noting that only terms higher than quadratic in 
R contribute after the differentiation, we obtain 



Hijki{x,x') 



1 



AM' 



{R^ + 3At2) ln(i?2 _ At 



2\2 



+ (^RAt 



~R~ 



In 



'R + At' 
R-At 



(A24) 



For convenience, we give the explicit forms for G'^xxx ^i^nd Ci^l^ here: 

^xxxx^'^J ) ^ n*^' ) "^^xxi^X 1 X ) -\- Hxxxx{x^ X ) 

1 



where 



127r2i?8a2 



(Ax' - Ar)(16Ax^ - 24Ax^At') - "iAt'R^ 



-(9At^ + mAx^Ae + lQAx^)Rt + (-72Aa;2At^ + 32Ax^At=^ + 32Aa;^)i?,^ 



At 



In 



R + At 



167r2/?9 \R-At 
-2AAt^Ax^Rl - SAx'^At^ + SAx 



R^t - (SAf + 9Ax')i?j^ 



R^t = Ay2 + Az' 
Ax = X — x' 

R^ = Rl + Ax^ = Ax' + Ay' + A^' 

a'^ = R^ - At' = Ax' + Ay' + Az' - At' 



(A25) 



(A26) 
(A27) 
(A28) 
(A29) 
(A30) 



To get G'^^}^^{x, x'), all we need to do is to replace i?' in Eq (|A25|) by R^ = Ay' + A 



X 
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